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ABSTRACT 

A system for the analysis of geomagnetic and gradient 
noise is described. The nominal 176 kHz Larmor frequency 
output of an optically pumped cesium magnetometer is mixed 
to a frequency in the 1.2 kHz range and FM tape recorded. 

The recorded output is applied to a periodmeter, the digital 
output of which is input to the XDS-9300 computer for 
analysis and display. The absolute accuracy of the period 
measurement is .75 ysec with a resolution of .25 ysec. 

The RMS noise of the tape recorder was determined to be 
.667 ysec or in the operating frequency range less than 
1 Hz or .28 gamma. Externally introduced field disturbances 
as small as .5 gamma were recorded and displayed, and the 
response time of the magnetometer determined to be less 
than 2.5 msec. 
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I. INTRODUCTION 



With a continuing interest in magnetic anomaly detec- 
tion the Naval Postgraduate School has studied geomagnetic 
noise, first with the proton precession magnetometer and 
later with the optically pumped rubidium type. In 1972 
the school obtained an optically pumped cesium magnetometer 
with accessories to convert it to a gradiometer. In order 
to take full advantage of the inherent high frequency 
response of this magnetometer a data collection system 
was developed to observe geomagnetic and gradient noise. 
Figure 1.1 is a block diagram of the system. 

A brief introduction to the earth's magnetic field 
and other measuring instruments is presented first followed 
by a description of the system developed. 
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Figure 1.1 System Block Diagram 
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II. THE MEASUREMENT FIELD 



A. MAIN FIELD 

World magnetic surveys both at ground level and more 
recently by satellite indicate that the principal source 
of the main geomagnetic field is beneath the earth’s crust. 
Approximately 90 % of the main field exhibits characteristics 
roughly equivalent to the field of a terrestrially centered 
short bar magnet or dipole inclined at 11.5° to the axis of 
rotation. See Figure 2.1. 




Figure 2.1 Dipole Appearance of Geomagnetic Field 
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The field strength is approximately 60 kilogaramas in 
the polar regions, 30 kilogammas in the equatorial region, 
and 50 kilogammas in the neighborhood of 37° North latitude 
on the west coast of the United States. Figure 2.2 shows 
a world chart of the total intensity of the surface field. ^ 
The main field has a slow secular change of about .1% per 
year and it is a regional phenomenon which may reach . 3 % 
in certain regions . 

B. UNITS 

The customary unit of measurement is the gamma, regarded 
interchangeably as a unit of induction (gauss) or of magnetic 
intensity (oersted) . 

1 gamma is 10”^ gauss (10"^ tesla) 
or 

— S —9 

1 gamma is 10”^ oersted (10“ /4ir ampere tum/meter) 

C. ELEMENTS 

The various elements of the magnetic field are shown 
in Figure 2.3. They are defined as: 

X North South Component 

^Air Force Cambridge Research Laboratories, Handbook of 
Geophysics and Space Environments , p. 11-11, 1965. 
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Figure 2.2 World Magnetic Chart 



Y East West Component 
Z Vertical Intensity 
D Declination 

H Intensity of the Horizontal Component 
F Total Intensity 
I Inclination or Dip 

ZEIIITH (-) 

.. t 

^ NORTH (+) 



Figure 2.3 Elements of the Field 




In terms of the more familiar bar magnet: A magnet 

perfectly free to turn in any direction would show the 
vector direction of the field. It is almost impossible 
to suspend a magnet like that so two are used. One with 
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a vertical axis (compass) and one with a horizontal axis 
(dip needle) . The compass shows the direction of the 
horizontal component H, and when properly oriented the dip 
needle shows the inclination or dip. 

D. TIME VARIATIONS 

The time variations of the geomagnetic field can be 
broadly categorized into quiet variation fields, disturbed 
variation fields, and micropulsations . 

1 . Quiet Variation Fields 

Geomagnetic records occasionally show a very smooth 
trace that clearly indicates a pattern of daily variation 
with respect to solar local time. This pattern is named 
"solar quiet day variation" and denoted Sq. There is also 
a lunar daily variation, denoted L, but it is much weaker 
and more variable. Both the Sq and L variations result 
from the d 5 mamo action of the atmosphere, a complex phenom- 
enon of convective motions caused by the solar heating of 
the upper atmosphere under the influence of the earth's 
magnetic field. These convective motions induce ionospheric 
currents and electric currents in the earth itself, both 
of which cause magnetic variations. The effects of electric 
current induction within the earth are sometimes referred 
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to as geological disturbances because the amount of induc- 
tion is dependent on the composition of the earth's crust. 
For example, the intensity of induced currents in the ocean 
is higher than those in the continental crustal layers. 

Figure 2.4 shows the Sq variation recorded in 
Monterey during a quiet period. 

2 . Disturbed Variation Fields 

Any fluctuation other than the quiet day solar and 
lunar variations can be called a magnetic disturbance. 

Other than disturbances directly attributable to a solar 
flare or fluctuations produced by upper atmospheric irreg- 
ular motions, all are caused originally by disturbed solar 
plasma and may therefore be generally classified as one 
family. The term magnetic storm is reserved to apply to 
a relatively severe, long lasting disturbance with recog- 
nizable features, but all other irregular variations 
deviating from the Sq and L variations are termed simply 
magnetic disturbances. They have durations lasting from 
100 hours to several minutes. Figure 2.5 shows a magnetic 
disturbance recorded in Monterey that immediately followed 
the Sq variations shovm in Figure 2.4. 
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Figure 2.5 Magnetic Disturbances, 100 Gammas Full Scale 



3 . Geomagnetic Micropulsations 



Geomagnetic micropulsations are field variations 
with periods ranging from .2 seconds to approximately 10 
minutes. A terminology recommended by the International 
Association of Geomagnetism and Aeronomy is pc for pulsa- 
tion continuous and pi for pulsation irregular. The period 
range for these are: 



pc 


1 


.2-5 sec 


pc 


2 


5-10 sec 


pc 


3 


10-45 sec 


pc 


4 


45-150 sec 


pc 


5 


150-600 sec 


pi 


1 


1-40 sec 


pi 


2 


40-150 sec 



A general picture of the geomagnetic field in the 
lower frequencies is shown in Figure 2.6. 

Pc 1 pulsations have very small amplitude, between 
.01 and .1 gamma and have a sinusoidal and beating form. 
The pulsations appear in bursts and usually last approxi- 
mately 30 minutes. They have been identified by their 
appearance and because there is a direct relationship 
between their reoccurrence period and mid- frequency . No 
seasonal or latitude dependence has been identified but 



21 




PERIOD^ (SEC) 



Figure 2.6 Geomagnetic Micropulsation Spectrum 

there is a diurnal variation in the mid-frequency events , 
high frequency events occur in post midnight hours, and 
low frequency events in the early afternoon. Figure 2.7 
is a record of a pc 1 micropulsation observed in Monterey. 

Figure 2.8 shows the typical form of the pc 2,3 
pulsations. They also exhibit a sinusoidal and beating 
character similar to the pc 1 t 5 rpe but their amplitude is 
slightly greater, usually .5 gamma. The duration of the 
signals seems to be related only to the sensitivity of the 
detector, which has lead to the h}rpo thesis that when no 
signal is observed the pulsation is not completely absent 
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Figure 2.7 Pci Micropulsation; 

1.8 gammas full scale; 45 second record 
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but rather too weak to be detected. Pc 2,3 events occur 
mainly during daylight hours and have a maximum near noon. 
The average period undergoes a diurnal variation, being 
a maximum at midday. They also have a 27 day (solar rota- 
tion) periodicity and a slight seasonal variation. 

As seen in the typical pc 4 record shown in Figure 
2.9, these micropulsations slowly rise to a maximum of 
between 5 and 20 gammas and then s}nmnetrically decay. The 
average signal lasts about one hour. In mid-latitudes a 
daytime occurrence pattern can be observed but the period 
is shorter in the daytime than at night, the opposite of 
the pc 2,3. 

The longest period variations, pc 5, display ampli- 
tudes from 10 to hundreds of gammas . Their form is again 
sinusoidal and they may last from several minutes to several 
hours. The largest amplitudes are observed in the auroral 
zone and the occurrences have a direct relationship with 
auroral zone phenomena. 

The class of micropulsations in the 2 to .005 Hz 
range with irregular form and a close relationship to upper 
atmospheric and auroral zone phenomena have been designated 
pi 1,2. Their amplitudes vary greatly, the larger sometimes 
being called micropulsation storms. Figure 2.10 shows the 
irregular appearance. 
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Figure 2.9 Pc4 Micropulsation; 

7 gammas full scale; 6 minute record 




Figure 2.10 Pi2 Micropulsation; 

7 gammas full scale; 6 minute record 
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III. MEASURING INSTRUMENTS 



Today a wide range of devices are available for measur- 
ing the magnetic field, ranging from the extremely old 
magnetized needle to the newer superconducting quantum 
types. The major considerations in the selection of a 
device in a non-observatory application are sensitivity, 
accuracy, frequency response, mobility, and reliability. 
Three types of magnetometers that have developed with these 
considerations as guidelines are the fluxgate, optically 
pvimped resonance, and superconducting quantvm. They will 
be described and a comparison of the different varieties 
and other t}q)es is given in Table I. 

A. FLUXGATE 

A fluxgate magnetometer, sometimes referred to as a 
saturable core or saturable inductor, measures magnetic 
fields by utilizing the variation of permeability with 
intensity of magnetization in a ferromagnetic core. An 
easily saturable ferromagnetic core of high permeability 
is driven cyclically into saturation by an applied magnetic 
field. In the absence of an external field, which is 
usually dc or very low ac and small with respect to the 
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applied field, the induced voltage in the sense windings 
will be syiimetrical, i.e., contain only odd harmonics of 
the fundamental driving current that produced the applied 
field. In the presence of an external field, the induced 
voltage becomes asymmetrical, the as3mimetry and the corre- 
sponding even harmonics being proportional to the amplitude 
of the component of the external field parallel to the axis 
of the core. 

The fluxgate magnetometer has been used widely in 
various applications where sensitivity, stability, and 
absolute accuracy has been sacrificed for the features of 
simplicity, low cost, low power, and ruggedness. Since it 
is a directional device, measuring only the component 
parallel to the axis of the core, a triaxial arrangement 
of three cores is necessary for total field measurements. 

B. OPTICALLY PUMPED RESONANCE 

The optically pumped resonance magnetometers (chiefly 
Rubidiiam, Cesixam, and Helixam) measure the total intensity 
of the magnetic field by an extremely accurate detection of 
the Zeeman effect. Due to the Zeeman effect the energy 
levels associated with a valence electron become split into 
sublevels whose separations are dependent on the total 
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intensity of the magnetic field. The detection of this 
splitting is based on the principles of optical pumping 
and monitoring. Optical pumping is described in a tutorial 
manner in Reference 4 and for the Cesium atom in Appendix A. 

Figure 3.1 shows a simplified block diagram of a typical 
pumping apparatus. A light source producing an intense col- 
limated beam of light is first filtered and circularly 
polarized and then passed through a resonance cell containing 
the vapor to be pumped. As the light is absorbed and emitted 
by the sample atoms, optical pumping proceeds and results 
in a steady state non-equilibriiam population distribution 
over the magnetic states of one of the ground levels. In 
this pxamped condition the sample is prepared for resonance 



LAMP 



RESONANCE 

LENS FILTER CELL 



LENS 




Figure 3.1 Typical Piomping Apparatus 
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and is characterized by a net magnetization M = along 
the ambient magnetic field direction. Resonance transitions 
between the Zeeman states of the pumped level are induced 
by subjecting the sample to an RF field of frequency 
sufficiently close to the resonant frequency fQ and with 
component perpendicular to the ambient field Hq = H^. 

The RF drive field introduces Zeeman coherence, and accord- 
ing to Bloch's equations the steady state behavior of the 
magnetization M changes and is no longer parallel to Hq but 
now has a transverse component Mj^ precessing around Hq at 
the frequency of H^. Consequently the component along 
Hq decreases, and the light transmission also. This varia- 
tion of light intensity with the applied RF, being a maximum 
at resonance, can be used as a detection method, or the 
photocell can be made to observe the signal corresponding 
to the M component, which is precessing or rotating at the 
same RF applied which induced the resonance. The photocell 
output is then at the resonant frequency. Appendix B 
details the two detection systems commonly used. The 
resonant frequency is proportional through a constant to 
the magnetic field (4.66 Hz/gamma for Rb; 3.499 Hz/gamma 
for Cs; 28 Hz/gamma for He) and thus the field measurement. 



29 



Figures 2.4 through 2.10 were all recorded using an 
optically pumped magnetometer and the appearance of the 
micropulsations differ somewhat from those one would find 
in a reference handbook. This illustrates an undesirable 
feature of this magnetometer when used in micropulsation 
work. It is equally sensitive to all frequencies and the 
larger amplitude, lower frequency changes such as the di- 
urnal variation continually drive the highly sensitive 
instrument off scale. The result is a sloping and discon- 
tinuous output that is difficult to interpret. 

C . CRYOGENIC SQUID 

A cryogenic SQUID (Superconducting Quantxam Interference 
Device) magnetometer consists of three basic parts: 1) a 

superconducting ring embodying a Josephson junction; 2) an 
externally generated HF field used to provide a small field 
variation superimposed on the ambient dc field and therefore 
to control the amount of flux through the ring; 3) an elec- 
tronic detection system to observe the impedance of the 
coil and make the output response a linear function of the 
ambient field. 

The operation utilizes the principle that if the flux 
through the ring is varied on either side of a minimum or 
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maximijm energy value, associated with the Josephson coupling 
energy, the SQUID will generate a field which will respec- 
tively oppose or add to the variation. The energy contribu- 
tion of the SQUID aiding or opposing the applied field 
variation according to whether the ambient dc field amounts 
to a minimum or maximum energy value constitutes a para- 
metric control of the inductance measured at the coil 
terminals. In other words, the SQUID alternately behaves 
diamagnetically and paramagnetically and will decrease and 
increase the inductive reactance. 

Through the use of a detector to monitor the coil imped- 
ance one can observe the variation of impedance with ambient 
magnetic field. Because the impedance variation is periodic 
a phase locked feedback loop is normally included to obtain 
a response linearly related to the field. 

Because SQUID compares the energy values of the ambient 
and weak superimposed fields, it measures only a change in 
field, not the total intensity. 
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D. GRADIOMETER 



A gradiometer , in theory, would measure the space 
gradients of the magnetic field, that is: 

^ ^ , dx 

3x 3y , 3z 

/ 3X > 9y 

3x 3y 3z 

^ , 3Z 

3x 3y 3z 



Only five of these are independent because X, Y, and Z 
are the derivatives of a potential and because of Laplace's 
equation. In practice gradiometers are constructed by 
taking the finite difference of two magnetometers, there- 
fore cancelling the ambient field that is approximately 
equal at both sensors. A local change or anomaly that is 
closer to one of the two sensors is not cancelled and 
therefore recorded. 

The disadvantage of the gradiometer is that with respect 
to a single sensor magnetometer most magnetic anomalies fall 
off as the inverse cube, 1/r where r is the distance from 
the anomaly to the sensor, but the gradiometer anomaly falls 
off as the inverse fourth power, l/r“^ where r is the 
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distance from the anomaly to the center of the sensor 
configuration. 

In Reference 10, Langan, L 3 mch, and Slack reported a 
gradiometer used in geophysical prospecting that measured 
the vertical gradient A H/ az where H was the magnitude of 
the anomaly and z the distance to the anomaly. It was 
constructed of two optically pumped sensors vertically 
separated and towed by a helicopter. 

Both the optically pumped and SQUID magnetometers are 
easily converted to gradiometers because of the form of 
their output, but they will measure different gradients. 
The optically pumped gradiometer, constructed of total 
field magnetometer sensors will measure the total field 
gradient AF/Az. A SQUID gradiometer, constructed of 
devices sensitive only to a field change will measure a 
field component gradient A Z/ Az, for example. 
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Table 1 Comparison of Magnetometers 



IV. SYSTEM DESCRIPTION 



The overall system can be broken into three subsystems: 

1) the data collection subsystem; 2) the periodmeter and 
interface subsystem; 3) the processing and display subsystem. 



A. DATA COLLECTION SUBSYSTEM 

The magnetometer used was the Varian V-4938 Cesium 
system. It used an optically pumped cesium sensor of the 
auto-oscillating type. A block diagram of the Varian sys- 
tem is shown in Figure 4.1. 



SENSOR 




Figure 4.1 Varian Magnetometer 
The sensor head produces the RF signal (Larmor frequency) 
directly proportional to the total magnetic field intensity. 
This signal together with the input power for the sensor is 
conducted along a 200 foot coaxial cable to the main readout 
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unit. The readout unit converts the Larmor frequency to a 
highly stable dc voltage which drives the pen recorder. 

The system was located in a wooded area approximately one 
mile from the Naval Postgraduate School, which provided a 
near noise free environment. 

A more complete block diagram of the V-4938 system is 
shown in Figure 4.2. The RF Larmor frequency, which was 
approximately 176 kHz, from the sensor is mixed in the read- 
out xmit to the 100 Hz to 3 kHz range and then applied to 
a highly stable frequency to voltage converter. This dc 
voltage, which now represents the field value, is passed 
through biasing circuitry to interface it to the pen recorder. 

The system configuration used for data collection differed 
from the Varian recommended system of Figure 4.2. This was 
necessary because of the requirement to record the output 
on magnetic tape. The actual system used is shown in Figure 
4.3. This configuration was used because the mixer output 
of the readout unit was a square wave and its frequency 
content exceeded the 5 KHz range of the tape recorder. The 
coupler unit shown connects between the readout unit and 
the sensor to provide dc power to the sensor and it also 
contains a mixer, but its output is a sinusoid rather than 
a square wave and recorded nicely. The coupler units primary 
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ELECTRONICS • H 



SENSOR HEAD 




Figure 4.2 Magnetometer Signal Flow 
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Figure 4.3 Magnetometer Data Collection Configuration 



purpose is to convert the system to a gradiometer which is 
shown in Figure 4.4. 




Figure 4.4 Gradiometer Data Collection Configuration 

Again this is not the configuration recommended but 
this provided both a magnetometer and gradiometer output 
suitable for recording and still allowed the pen recorder 
to monitor the gradiometer signal. The bias applied to 
the one sensor was required to generate a difference 
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frequency between the sensors sufficient to drive the mixer 
circuitry. The sensors were located approximately 180 ft 
apart resulting in a Larmor frequency difference of only 
approximately 30 Hz, and the mixers required a difference 
of at least 100 Hz and optimally between 1 and 3 kHz. 

Figure 4.5 shows the readout units mixer output with an 
insufficient difference frequency and Figure 4.6 shows the 
coupler units sinusoidal mixer output. 

One other system characteristic imposed a limitation, 
that was the orientation sensitivity of the single cell 
sensor. For optimxam performance the sensor head needed to 
be oriented both toward magnetic North and 45° to the total 
field vector. Figure 4.7 shows the non-metallic housing 
and Figure 4.8 the geometry. 

Figure 4.9 shows the system components and equipment 
used to monitor and record the data. 
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Figure 4.5 Distorted Readout Unit Mixer Output 




Figure 4.6 Coupler Unit Mixer Output 
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Figure 4.7 



Non-Metallic Housing Used to Fulfill 
Orientation Requirements 
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field 




Figure 4.8 Sensor Housing Geometry 
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B. PERIODMETER AND INTERFACE SUBSYSTEM 



The frequency response of the magnetometer when used 
with the readout unit is limited to 20 Hz by the discrimin- 
ator and even more by the pen recorder response. The 
theoretical frequency response of the sensor is limited 
only by the bandwidths of the amplifier and feedback loop, 
which are extremely wide. To take advantage of this the 
mixer output was recorded on an HP -3960 FM tape recorder. 

This signal still contained the high frequency components 
present in the Larmor signal. A periodmeter and interface 
unit were constructed to measure the period of successive 
cycles and interface this digital count into the XDS-9300 
computer, for further display on the ACT -10 graphics console. 

Figure 4.10 is a block diagram of the periodmeter. The 
data signal is passed through a simple RC high pass filter 
to reject any 60 Hz modulation that might be present and 
then input to a zero crossing detector. The falling edge 
of this is used to trigger a one-shot that inhibits the 
otherwise free running 4 MHz count. The one-shot pulse 
is inverted and delayed and the rising edge used to trigger 
flip flop type memories. Once the count is latched another 
one shot is used to clear the count and provide a signal 
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Fig 4.10 PERIODMETER 



indicating the cycles completion. The cycle then repeats. 

The actual circuits and timing diagrams are presented in 
Appendix C. 

A direct access channel into the XDS-9300 was used to 
input the binary count directly into sequential memory loca- 
tions. Figure 4.11 is a block diagram of the periodmeter 
and interface. To implement the input the following sequence 
of instructions were used: 

EOM 032400 
SKS 032400 
BRU $-1 
PIN ITEM 

With the execution of the EOM instruction the contents of 
the control register are available at the output connector. 
This is decoded and toggles a flip flop to alert the period- 
meter that the XDS is ready to accept an input. This EOM 
alert and the periodmeter ' s completion signal are logically 
combined and send a ready signal to the XDS-9300. Upon 
receipt of this ready signal and the execution of the PIN 
the data is transferred and stored in the location specified 
in the PIN instruction. The SKS instruction is used to sense 
if both the computer and periodmeter are ready. If both the 
EOM alert and completion signal are not present, the next 
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Figure 4.11 Periodmeter and Interface Block Diagram 



s . 



instruction, branch unconditionally back to the SKS , is 
executed. This looping continues until both systems are 
ready. If both are ready the branch instruction is skipped. 
The ready signal is also used to gate the data through in- 
verters because it is inverted prior to being transferred 
into the XDS-9300. The actual circuit diagrams are pre- 
sented in Appendix C. 



C. DISPLAY SUBSYSTEM 

A block diagram of the data flow in the interface and 
display subsystems is shown in Figure 4.12. 




ACT- 10 

GRAPHICS 

CONSOLE 



Figure 4.12 Data Flow From Periodmeter Through Display 



The physical block length of data input to the XDS-9300 
was 12,288 points. A mean and standard deviation of these 
points was found and values not within five standard 
deviations were discarded. This eliminated erroneous data 
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points caused by noise spikei^ either in the computer, 
periodmeter, or on the recording. The block was then 
buffered out in groups of 4096 points and written on 
magnetic tape. To be displayed the data was buffered 
back in, again in groups of 4096 points, normalized, and 
displayed on the AGT-10 graphics console. Both programs 
are self-explanatory and are presented in Appendix D. 
Signal flow diagrams for the two are shown in Figure 4.13 
and 4.14. 
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Figure 4.13 Signal Flow Diagram of Input Program 
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Figure 4.14 Signal Flow Diagram of Display Program 
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V. RESULTS AND CONCLUSIONS 



The system was tested in three phases: 1) laboratory 

signals; 2) tape recorded laboratory signals; 3) magnetom- 
eter signals. 

First, constant frequency signals from a laboratory 
generator were applied. The frequencies sampled and the 
resulting period measurement are listed in Table 2. 

FREQUENCY Hz MEAN PERIOD ptsec THEORETICAL PERIOD ^sec 



1000.0 


999.25 


1000.00 


1000.1 


999.00 


999.00 


1000.5 


998.75 


999.50 


1001.0 


998.25 


999.00 


1003.0 


996.25 


997.01 


1005.0 


994.25 


995.02 



Table 2 



The .75 p sec difference between theoretical and actual is 
due to the clear and latch time of the periodmeter. These 
measurements were within +1 count or +25 psec with a 4MHz 
time base. This established the resolution of the device 
as .1 Hz or better. The outputs displayed for the sampled 
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frequencies of 1000.0 Hz, 1000.1 Hz, and 1001.0 Hz are 
shown in Figures 5.1 - 5.3 and a distinguishable difference 
is noticeable. To insure the system would follow a rapidly 
changing frequency a swept frequency signal was applied and 
the display shown in Figure 5.4. 

The second phase of testing was designed to determine 
the magnitude of the tape recorder flutter and noise. It 
involved recording the same frequencies used in phase 1 and 
comparing the results. The frequencies, mean period measure- 
ment, and standard deviation for each are listed in Table 3. 

FREQUENCY Hz MEAN PERIOD usec STANDARD DEVIATION ^sec 



1000.0 


999.25 


.575 


1000.1 


999.00 


.700 


1000.5 


998.75 


.850 


1001.0 


998.25 


.550 


1003.0 


996.25 


.650 


1005.0 


994.25 


.675 



The mean period, calculated for 12,288 periods, remained 
the same. The average standard deviation was .667 psec. 
This established the RMS noise level of the system when 
operated in the 1200 Hz range to be approximately 1 Hz or 
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Figure 5.1 Output Display of Direct Input 
1000.0 Hz Laboratory Signal 



1000. 



PERIOD 

pSEC 




Figure 5.2 Output Display of Direct Input 
1000.1 Hz Laboratory Signal 
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1000 




F igure 5 . 3 



Output Display of Direct Input 
1001.0 Hz Laboratory Signal 




TIME, SEC 



1000. 



PERIOD 

jUSEC 



990. 

0 . 



Figure 5.4 Output Display of Direct Swept Frequency 
Laboratory Signal 



56 



.27 gamma. The AGT-10 displays for the sampled frequencies 
of 1000.0 Hz and 1000.1 Hz are shoxm in Figures 5.5 and 5.6. 
The previously discernible difference is now buried in tape 
recorder noise. Figure 5.7 shows the frequency transition 
from 1001.0 to 1003.0 Hz and this 2 Hz difference is easily 
distinguishable. 

The final phase of testing involved introducing a known 
magnetic disturbance and analyzing the magnetometer output. 

A 72 turn Helmholtz coil, in the configuration of Figure 5.8, 
was used to introduce the disturbance. The sensor was 
placed on the axis of the coil and the coil aligned with 
the total field vector. This insured the induced field of 
the coil would directly aid or oppose the ambient field 
and allowed changes of predetermined amounts. With the 
simplified geometry 



H = 



N I R" 



4 IT 



10 






where N = number of turns in the coil 



I = current in the coil 



R = radius of the coil 



r = distance from coil to sensor 




was required for conversion from 
ampere-tums per meter to gammas. 
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990 




0 . 1.2 

TIME, SEC 

Figure 5.5 Output Display of Recorded 
1000.0 Hz Laboratory Signal 



1000. 
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^SEC 



990. 




Figure 5.6 



Output Display of Recorded 
1000.1 Hz Laboratory Signal 
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TIME, SEC 

Figure 5.7 Output Display of RecordedlOOl. 

and 1003.0 Hz Laboratory Signal 
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Figure 5.8 Coil and Sensor Geometry Used 
to Disturb the Field 
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J 




The coil utilized had a mean radius of 33 cm and was placed 
3.5 m from the sensor. Table 4 tabulates the various 
currents and observed versus predicted changes. Figures 
5.9 through 5.30 show the magnetometer pen recorder output 



and a 


selected time 


frame of the 


corresponding 


AGT-10 display 




THEORETICAL 


MAGNET. 


THEORETICAL 


DISPLAY 


I ma 


AH 


AH 


Af 


Af 


400 


45.964 


46.0 


160.874 


160.893 


300 


34.473 


35.0 


120.655 


120.666 


240 


27.578 


27.5 


96.524 


96.531 


200 


22.982 


22.7 


80.437 


80.442 


160 


18.386 


18.7 


64.350 


64.354 


120 


13.789 


13.8 


48.262 


48.264 


100 


11.491 


11.1 


40.218 


40.219 


90 


10.342 


9.6 


36.197 


36.198 


80 


9.193 


9.0 


32.175 


31.176 


70 


8.044 


8.4 


28.153 


28.154 


60 


6.895 


6.8 


24.131 


24.132 


40 


4.596 


4.2 


16.087 


16.086 


20 


2.298 


2.1 


8.044 


8.045 


10 


1.149 


1.1 


4.022 


4.023 


5 


.574 


.65 


2.011 


2.012 



Table 4 
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FIELD CHANGE 




Figure 5.9 Pen Recorder Output for AH= 45 gamma 




TIME, SEC 

Figure 5.10 AGT-10 Display for AH= 45 gamma 
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FIELD CHANGE 




Figure 5.11 Pen Recorder Output for AH= 35 gamma 
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Figure 5.12 
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TIME, SEC 



Figure 5.13 Pen Recorder Output for A H= 28 gamma 
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Figure 5.14 



AGT-10 Display for 



AH= 28 gamma 



64 



FIELD CHANGE 




Figure 5.15 Pen Recorder Output for aH= 18 gamma 
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Figure 5.16 



AGT-10 Display for 



A H= 18 gamma 
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Figure 5.17 Pen Recorder Output for aH= 11 gamma 
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Figure 5.18 AGT-10 Display for A 11 gamma 
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FIELD CHANGE 




TIME, SEC 



Figure 5.19 Pen Recorder Output for A H= 8 gamma 



1000 , 




0 . 1.2 
TIME, SEC 

Figure 5.20 AGT-10 Display for A H= 8 gamma 
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FIELD CHANGE 




Figure 5.21 Pen Recorder Output for A H= 4 gainma 
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Figure 5.22 AGT-10 Display for AH= 4 gainma 
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Figure 5.23 Pen Recorder Output for AH= 2 gamma 
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0 . 1.2 
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Figure 5.24 AGT-10 Display for AH= 2 gamma 
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FIELD CHANGE 




Figure 5.25 Pen Recorder Output for AH= 1 gamma 



1000. 




TIME, SEC 

Figure 5.26 AGT Display for AH= 1 gamma 
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FIELD CHANGE 




Figure 5.27 



Pen Recorder Output for AH= .6 gamma 



1000. 




0 . 1.2 

TIME, SEC 

Figure 5.28 AGT-10 Display for AH= .6 gamma 
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To obtain a precise measurement of the magnetometer's 
response time the AGT-10 display was expanded. Figure 5.29 
shows the rising edge and Figure 5.30 the falling edge. As 
seen the rising edge required approximately three periods 
or 2.5 msec to respond and the falling edge only one period 
or .8 msec. 




TIME, SEC 

Figure 5.29 Expanded View of Rising Edge 
of Magnetometer Response 
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TIME, SEC 

Figure 5.30 Expanded View of Falling Edge 
of Magnetometer Response 
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The periodmeter and interface unit measured the period 
of the signal with an absolute accuracy of .75 |jl sec and 
properly input the value to the computer. The resolution 
of the device was .25 sec using a 4 MHz time base and 
this could be increased by using a faster time base. The 
system was successfully operated with a 10 MHz clock fre- 
quency and resolution of .1 |j.sec. The use of an increased 
clock frequency or averaging of periods to increase the 
periodmeters resolution would not improve the system sensi- 
tivity. The overall system sensitivity is limited by the 
flutter and noise of the tape recorder. TV'/o methods were 
attempted to overcome this: 1) the use of a s}mchronized 

oscillator constructed of two multivibrators synchronized 
with the zero crossings of the input signal; 2) recording 
a 1 kHz reference signal on a second channel of the 
recorder and multiplying it up for use as a time base. 

The pulse width instability of the multivibrators was 
greater than the flutter noise so the synchronized oscillator 
was unsatisfactory. A phase lock loop t 5 q)e frequency multi- 
plier was constructed to multiply the 1 kHz to 1 MHz but 
the delay time through the loop was too great to cancel the 
flutter for a period to period measurement. 

To achieve greater system sensitivity and resolution some 
means of flutter cancellation is necessary. 
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APPENDIX A 



OPTICAL PIMPING OF THE CESIUM ATOM 



Figure A-1 shows the energy levels of the cesixam atom. 

Absorption of photons of light by the cesium valence elec- 

2 

trons cause their energy to be increased from the ground 

2 

state to the Pi state. Due to the light's polarization 

the magnitude of the quantxim number m is also increased by 

one. An excited atom quickly emits its energy and returns 

to any of the ground state sublevels. The absorption and 

re-emission is repeated resulting in a shift to the m=+4 

2 

sublevel. Since the Pi state contains no m level above 

2 

+4, transitions out of the Si , m=+4 level are not possible. 
This is precisely optical pumping and under these conditions 
there can be no further absorption of light. 



75 



ENERGY 

FREQUENCY 

m»*4 



m-o 



m— 4 
m--3 



n-o 



ma+3 




MAGNETIC FIELD STRENGTH 




Figv.re A-1 Energy Levels of Cesium Atom 
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APPENDIX B 



OPTICAL DETECTION SYSTEMS 



The first method used to obtain a continuous record of 
the total field intensity used a control loop arrangement 
and was therefore called a Controlled Magnetometer, A 
block diagram of it is shown in Figure B-1. 



DIRECTION OF FIELD 

► 



FILTER resonance 

CELL PHOTOCELL 




A small frequency modulation is applied to the local RF 
oscillator to sweep it back and forth through resonance. 
Owing to this sweeping the signal collected at the photo- 
cell is also modulated at the same frequency, but its phase 
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( 



with respect to the initial modulating signal depends on 
the position of the local RF with respect to the resonant 
frequency. The photocell output and modulating signal are 
phase detected and the output is used as an error signal 
in a control loop to maintain the local RF oscillator coinci- 
cent with the resonant frequency. This system does not place 
severe restrictions on the performance of the photocell or 
its amplifier, has high accuracy, good stability, and good 
low field capability. 

One attribute the controlled system does not have is a 
good response to transients. Figure B-2 shows the block 
diagram of a system that has the excellent response to 
transients necessary to detect and observe micropulsations. 



DIRECTION OF 
FIELD 




Figure B-2 Self-Oscillating Magnetometer 
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In theory a second light beam is placed at right angle with 
respect to the first and the variation of the M component 
of the magnetization is observed. Since this transverse 
component oscillates at the frequency of the signal injected 
and the amplitude is a maximum at resonance, the photocell 
is simply amplified, phase shifted, and re-injected to the 
coils. The system will then oscillate at the resonant 
frequency. In practice it is possible to use a single light 
beam inclined at 45° to the direction of the field to be 
measured. This places an orientation restriction on the 
sensor, but systems employing two and three resonance cells 
have been developed that overcome this. This system places 
severe restrictions on the amplifier and especially on the 
accuracy of the phase shift, which is required to maintain 
oscillation. It has a low absolute accuracy and poor long 
term stability but as previously mentioned, excellent 
response to transients. 
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APPENDIX C 



PERIODMETER AND INTERFACE CIRCUITS 
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Figure C l Counting Circuitry 
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Figure C-2 Input Circuitry Timing Diagram 
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Figure C3 Input Circuitry 
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FIGURE C4 TIME BASE 
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FIGURE C 5 COMPUTER INTERFACE 
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APPENDIX D COMPUTER PROGRAMS 



-JOE 

-FCF755/N LS,GC,S 

CI^'E^SICN CGUNT(13CC0) ,!VAR(6) 
IMcGER CCLNT 
CATA NSAMP, MULT/4096, 3/ 

LELK = NSAMP’»MULT 
106 FCPMAK IX, 12( 110) ) 

999 CCNTINLE 

INITIALIZE VARIAL8LSS 

INA)«=1 
IMN = ICOOOO 
SL^=0. 

SlMSC=0. 

C = C. 

IEAC=0 
EAC=C. 



INPUT CATA FROM PERICCMETER/ INTERFACE 



CC ^ l = i,LI 
ECM C-2400 
SKS 032400 
EFL $-1 
FIN ITEM 
ICCLNT =0 

ICCLNT=LANC( ITEM, 177B) 

ICCLNT=LICR(IC0UNT,LANC(LRS(ITEM,1) ,37600 6) ) 
ICCLNT=LICR(ICOUNT, LANC(LRS( ITEM ,1) ,1400008 ) ) 

CCLNTd ) = ICCUNT 

2 CCNTINLE 

FCFM AN INITIAL AND THEN PRECICE MEAN £ SDEV 

CC 3 I=1,L6LK 
SIM*SUM+CCLNT (I ) 

SLMSC=SLKSC: + COUNT (I )^=»2 
C=C4l.O 

3 CCNTINLE 
MEAN=SLM/C 

SCEV=SQRT IA8S(C*SU«SG-SLM=S'«2) )/C 
VPITE (6 ,100 )MEAN,SCEV 
ICC FCFMAT(»0 MEAN= ‘ , F 12. 4/ 

1* SCEV=',F12.4) 

CC 4 I=1,L8LK 

1F(ABS(CCUNT( I)-MEAN).LE.5.*SC£V) GC TO 4 
SLM=SUM-CCLNT (I ) 

StMSG = SLMSC-COUNT (I 
E AC = 8AC-H. 

4 CCNTINUE 
C=C-EAC 
MEAN=SUM/C 

SCEV=SQRT ( AES(C*SLMSC-SLM=«''-i'2) )/C 
VFITE(6,101)MEAN,SCEV 

101 FCFMATCO TRUE MEAN= » , F 12 .4/ 

1' TRUE SDcV=* ,F12.4 > 

FCFM A NEk ARRAY CONSISTING OF POINTS WITHIN 10 SDEV 

IEAC=0 

N£VI=0 

CC 6 I=1,LBLK 

IF (ABS(CCLNT(I )-MEAN).LE.10.=<'SDEV)NEk = NEkI+l;GG 7C 7 
I EAC=I8AC4l 

kRITE(6,lC2) I8AD, I , CCLNTd ) 

102 FCFMAK’O THIS IS BAC PCINT NO', 15/ 

1* ITS OLD SUBSCRIPT WAS', 18/ 

2' ITS VALUE WAS', 110) 

GC TC 6 

7 CCLNKNEWI )=COUNTd ) 

6 CCNTINUE 



I=1,L8LK 
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CC 15 I=NEU+1,LBLK 
CCCMd ) = C 
15 CCMINUE 

kFITE (6 ,103) NEW! 

10 - FCFNAK *C NEWI=* ,110) 

ZERO FILL 4SRR/1Y FROM LAST POINT TO ENC OF BLOCK 

CC 6 I=NEWI+1,L6LK 
CCLNTd )=0 
£ CENTIME 

FINC A MINIMLP AND MAXIMUN VALLE 

CC S I*1,NEWI 
ININ=AMINO<CCUNT( I ) ,IMINl 
IPAX=AMAX0(C0UNT( I), I MAX) 

S CCNTINLE 

WRITE(6 ,10A)IMIN,!NAX 
lOA FCFMATCC MIN IMUM= • , 1 1C / 

1* MAXIMUM^' ,110 ) 

ELFFEFCUT TG TAPE THE ARRAY IVAR CONTAINING VAFIA6LES 

UAP(1)=NEUI 
I\AR(2)=MEAN 
I VAR(3)=SCEV 
I\AR{4)=IMIN 
IVAF<5) =IMAX 

WFITE(6,105) ( iVARd ) ,1=1, 5) 

1C5 FCRMAT( IX, 5110) 

CLTFLTl lOi) 'WRITE ON TAPE* 

REAC( 1C1,601)IANSS 
6C1 FCFMAT(A4) 

IF(IANSS.NE.3HYES) GO TG 603 
CALL BLFFERCUT( 1, IdVAR ,5 ,INC) 

11 IFdNCoEQ.l) GO TO 11 

ELFFEPCUT TG TAPE THE ARRAY CGLNT IN GROUPS OF 4096 

CC 12 1=1, MULT 
J = I-1 

K = ( J^NSAMF ) + l 

CALL 6UFFERCUT(1,1,CCUNT(K) ,N£AMF,INC) 

1 ' I F ( INC. EC. 1) GO TO 13 

12 CCNTINLE 
6C- CCNTINLE 

CLTFUT(1C2) 'READ ANCTHER 1A336' 

F£AC( 1C1,600 )I ANS 
6CC FCFMAT1A4) 

IF dANS.EC.SFYES ) GO TO 999 
ENC 

^LGAC XF,MAP 
^CATA 
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OOr-k ooo 



IPAICI- 

0>C/T^ 

007067 

$PNC 

-^G7 

^JCE 

-nFCPTP4h 



0C1C6711 



LS,GC 

Cl^' ENSIGN I AX IS (10), INCAT A (4105) ,CATA(410' 
II CCIR(5 ), ITCIR( 15 ) , ITEXT(24) ,T(2 ) 

INTEGER RSSP(4) 

CATA NSAHP,MULTtNCIS/4CS6,3, 16/ 

IC£V=1 

LELK=NSAMF#^'L•LT 

CIV*FL0AT(MCIS) 

NELK=1 

lAXISt-l )=IhEAD(0, 10) 
I/XlS(2)=IPACK(-.5,.75,0) 

I4>IS(2 )=IPACK(-.5tl.,0 ) 

IA>IS(4 )=IPACK(-.5,0.,1) 

I4XIS(5 )=IFACK(0.,C., 1) 

I/XIS(6 )=IPACK(0. ,-.01,1) 

IAXIS( 7 )=IFACK(0.,0.,0J 
IAXIS(8 ) = IFACK( .5,0., 1 ) 
IAXIS(9)=iPACK(.5,-.Ql,l) 

999 CCNTINCE 
TINE=C. 

LFEC = NSA^^P 

2 CALL FNS( ICEV , I S W , I ER ) 

If { I6R.NE.G) 

1F( ISW.NE.Q) 



fionfitKJ 

CUTPtT(101)I5R,'Fl* 
GO TC 2 



BCFFEFIN ARRAV IVAR kITH VARIABLES 

CALL EU FFEPIN (1, 1, IVAR, 5, INC ) 

3 IF(INC.EC.l) GO TC 3 
NEkI=lVAR( 1) 

NE AN = FLCAT( IVAR(2 ) )/4. 
3CEV=FLCAT( IVAR(3 ) )/4. 
r'IN = FLQAT(IVAR(4) )/4« 
yAX=FLGAT(IVAR(5) )/4. 
kFITE(6,8Q0)( IVARd ) ,1=1,5) 

8CC FCFNATdX, *N0. CF PC I NT S= ' , 1 10/ 
1' N£AN=',I10/ 

2‘ SEEV=',I10/ 

3' MNIMLM=‘ , 110/ 

4* NAXIMLM=* ,110) 

kRITE NEAN ,SCEV ,MIN, yAX CN SCREEN 

CALL CGINITdOEV, IGCIR,5,IEP) 

IF( IER.NE.0)0UTPUT( 101 )IER, ‘CGI* 
CALL DTIN!T(IDEV,ITCIR,15,IER) 
IFdER.NE.C) OUTPLT(lOl) I5R,*D2 
PNrrnP(<;A. TTpyr j 



2C0 



1,3,IER) 

, 



irllcK«TNu«Ul UUirLIl 

ENCCD£(96,200,IT£XT) 

FCFNAT( 'PERIOD’ ) 

CALL TEXTC(IDEV,ITcXT,24,3,26, 1,3, 

IF( 1ER.NE.0)0L'TPLT( 101) lER, 'PER ICC 
ENCCCE( 96,201, ITEXT) 

201 FCF^AT( *T SEC ) 

CALL TcXTG(IDEV,IT5XT,24,22,75,l,3,IER) 
IF ( 1 ER. NE.C ) OUTPUT! 10 DIER, 'PTINE' 



30: 



3C4 



306 



CALL TcXTGdDEV, ITEXi ,24,2i 
IF ( 1 ER. NE.C ) OUTPUT! 10 DIER. 

ENCCDE! 96 ,303, ITEXT)NBLK 
FCFNAT! 'STATISTICS ON BLOCK', 

CALL T S XTC( IDEV, ITEXT, 24, 24,3 _ 

IF (lER.NE.O) OUTPLTdODIER, 'TEl* 
£NCCDE(96,3C4, ITEXT)MEAN,SDEV 
FCFNAT! 'MEAN VALUE= * , F8 .3 , 4X , ' S C 6 V = » , F8 
CALL TE XTO d DE V, I T6 XT ,24,2 5,24,1 ,2 ,! ER) 
IF! lER.NE.O) OUTPUTdOl )IER, 'T82' 



,13) 

32, 1,3,IER) 



F8.3 ) 



IFdER.NE.O) OUTPUTdOl 

ENCCC£(96, 206, ITEXT) NIN,MAX 
FCFMAT!'MIN PERIOD IN MICSGSEC 
14X,*MAX PERIOD=' ,F8.3 ) 



- I 



F8.3 



S7 



,IVAR !6 ) , 




I 

I 

1 



I 

■! 




orio ooo ooo ooo 



208 
2 10 
209 

^0 

41 

12 C 

215 

218 

998 



C 4 LL 7 E>TC(IDEV,ITEX 7 , 24 » 26 , 24, 1 , 2 ,ieS) 

IF (lER.NE.O) OU 7 PUT( 101 )IER, ' 7 E 2 ‘ 
E^CCDE(96,2C8,!7EXT) 

FCF^A 7 ('S 67 MIN £ MAX VALUE CF Y AXIS IN MICRSEC) 



24 , 1 , 
• 7 E 4 * 



2 ,IER) 



CALL 7 SXT 0 (ICEV, 178 X 7 , 24 , 28 , 

IF(IER.NE.O) OU 7 FUT(lOl )I£R, 

ENCCDE( 96 , 210 ,I 7 EXT) 

FCFMA 7 (' 7 YPE A MIN, A COMMA, AND MAX IN AN F FCPMA 7 M 
CALL 7 £X 7 Q(IDEV,I 7 EXT, 24 , 29 . 24 , 1 , 2 ,IERI 
IF( 1 ER.NE. 0 )GUTPU 7 ( 101 )IER, ‘PTYPE ' 
ENCCCE( 96 , 209 ,I 7 EX 7 ) 

FCRMA 7 ('TQ USE STA 7 IS 7 ICAL MIN ANC MAX AS Y AXIS 
IMIN ANC MAX TYPE JUS 7 A C/R') 

CALL ! 5 > 70 nDEV,I 7 FXT, 24 , 22 , 24 ,l, 2 ,IER) 



IF(IER^NE.C) OU 7 PL 7 ( 101 ) IER,*T 1 



24 , 

I 



1 , 2 ,IER ) 



CALL 7 EX 7 R(I 06 V, 77 777 7778 , 1 , 20, 2 
IF( IER.NE. 0 )QU 7 PU 7 ( 101 )IER, 'TXTR 
1 F(MCC(ITD!R( 9 ) , 8 ).EQ.O) GO 70 40 
CALL 7 EXTI(IDEV,RBSP, 4 , 0 , 9 , lER) 

I F (PESP m .N 6 . 0 H > 6 C 70 41 
CMN = MIN 
CM AX=MA> 

CCN 7 INLE 

CECCCEI 16 , 120 ,R£SF )CMIN,DMAX 
FCFMA 7 ( 2 F 8 . 0 > 

1 F(CMIN .GE.DMAX) GC TO 42 
CM AX 2 =CMAX-CMIN 
ENCGCEI 96 , 215 ,I 7 EX 7 )DMIN 
FCFMA 7 ( F 8 . 2 > 

CALL 7 EX 7 C(ICEV, I 7 £XT, 24 , 20 , 20 , 1 , 2 ,IER) 
IF( IER.NE.C) 0 UTPU 7 ( 101 )IER, • LAS 7 • 
EMCCCE( 96 , 216 ,I 7 EX 7 )DMAX 
FCFMA 7 (F 8 . 2 ) 

CALL 7 EX 70 (IDEV,I 7 EX 7 , 24 , 5 , 

IF (IER.NE.O)OUTPU 7 < 101 ) lER, 

KKK = 1 
CCMINLE 



20 , 1 , 2 , lER) 
• LAST + 1 » 



ELFFERIN 4096 PQIN 7 S CF ARRAY CCUN 7 IN 7 G INCA 7 A 
IC£V =1 

CALL eUFFERIN(l,l,INDA 7 A,NSAMP, INC) 

4 IF ( INC. 8 C. 1 ) GO 7 G 4 

I F ( (N 8 LK.LE.il ) .OR. (NELK.GE. 14 )J GC TO 301 
CC 27 Ll=l, 4081,16 

WFI 7 E( 6 , 106 M INOAT A ( 1 2 ) , L 2 *L 1 , L 1+ 15 ) 

1 C 6 FCFMAT( IX, 16 ( 16 , 2 X ) ) 

27 CONTINUE 
8 C 1 CCM 7 INUE 

MAKE INCA 7 A A REAL ARRAY 

CC 10 I= 1 ,NSAMP 
CA 7 A(I)=FLCAT(IN 0 ATA(I ) )/ 4 . 

1 C CONTINUE 

IF ( (KKK.EC .3 ) .AND. (NEV^I.NE.LBLK ) ) LR £C=N EW I - 2 ’tNS AM F 
NORMALIZE (0 TO 1 ) TC CISPLAY 
CC 12 I= 1 ,LREC 

C ATA (I ) = (OATA( I )-CMIN ) /CMAX 2 
12 CONTINUE 

FCFM TME AFRAY INOATA GRAPHICS V^CRDS 



L = 1 

996 CCNTINLE 

CALL CGINlTdDEV, IGCIR, 5 ,ISR) 

IF(IER.NE.O) OUTPUTdOl) IER,‘D 1 ' 

7 (1 )=TIME 
CC 19 !» 1 ,NDIS 

T IM E = T I ME+( { (CATA (L)*CMAX2>+DM IN )/ 1000000 . ) 
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ooo 



IF(I.EQ.(NCIS/2) >T(2)=TIME 
1F< I.EQ.NDIS) TO)=7IME 
15 CCMINLE 

ENCCCE( 56 ,320,I7EXT )T(1>,T(2J,T(2) 

320 FCFMAT(28X,F6,3, 12X,F6.3,12>tF6.2) 

CALL TEXTOilDEV, ITEXT,24,22 ,1,1»3,IER) 

IF( lER.NE.OOUTPUT (101) lER, 'TTIVE* 

UCATA( 1) = II-EAD(0, 10) 
INCATA(2)=IPACK(-.5,.5,0) 

I^CATA( 3 ) = IPACK( ( 1 . /D I V )-. 5 ,CAT A ( L ) ,0 ) 

CC 6 I=2,NCIS 

INCATA( 1+2) = I PACK ( ( I /D I V )-. 5 , DATA ( L + I-1 ) ,1 ) 
6 CCMINUE 
NVCFD=NCIS+2 

CFAk CN SCREEN 

CALL GRAPHCCIDEV, lAXIS.S,!, lER) 

IF( 1ER.NE.O)OUTPUT(101)IER, ’Gl» 

CALL GRAFHCdCEV, IN C AT A .NWORC , 2 , I ER ) 
IFdER.NE.OQUTPUTdODIER, 'G2* 

N=N+NCIS 

eC CALL FNS ( ICEVtISw ,IER) 

IF(IER.NE.O) OUTPliTdOl) IER»'F1' 

IF (LLS( ISU,15).LT.O ) GC TG 18 
IF(LLS( ISW,16).LT.O) GG TC 995 
C-C TO 80 
1€ L = L + NL'IS 

IF{L.LE.(NSAMP-NDIS+1) ) GO TC 996 
Kt<K = KKK + l 

IF(KKK.C-T,3 )N8LK = NELK + 1 ;G0 TO 999 
GC TO 998 
995 STCF 
£^C 

*LCAC XF,KAP 
>CATA 
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